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Ultrahigh energy neutrinos (UHEz/) scatter on relic neutrinos (R^) producing Z bosons, which can 
decay hadronically producing protons (Z-burst). We compare the predicted proton spectrum with 
the observed ultrahigh energy cosmic ray (UHECR) spectrum and determine the mass of the heaviest 
Kv via a maximum likelihood analysis. Our prediction depends on the origin of the power-like part 
of the UHECR spectrum: m v = 2.75±J;;j? eV for Galactic halo and 0.26t£?° eV for extragalaetic 
(EG) origin. The necessary UHEj/ flux should be detected in the near future. 

PACS numbers: 14.60.Pq, 98.70.Sa, 95.85.Ry, 95.35.x+d 



I. Introduction. — The interaction of protons (p) with 
photons (7) of the cosmic microwave background radia- 
tion (CMBR) predicts a sharp drop in the cosmic ray 
flux above the Greisen-Zatsepin-Kuzmin (GZK) cutoff 
around 4 • 10 19 eV |Q. The available data show no such 
drop. About 20 events above 10 20 eV were observed by 
experiments such as AGASA Fly's Eye j|, Haverah 
Park |J, Yakutsk §], and HiRes g. The attenuation 
length of protons above the GZK cutoff is « 50 Mpc; 
but no obvious astrophysical source candidate is known 
within this distance. No conventional explanation for the 
observed UHECR spectrum is known (7). 

Already in the early 80 's there were discussions that 
the UHEj/ spectrum could have absorption dips at ener- 
gies around E T v f = M|/(2m„J = 4.2 ■ 10 21 (1 eV/m v .) 
eV due to resonant annihilation with the R^s, predicted 
by the hot Big Bang cosmology, into Z bosons of mass 
Mz Recently it was realized that the same anni- 

hilation mechanism gives a possible solution to the GZK 
problem |l(J. ft was argued that the UHECRs above the 
GZK cutoff are mainly from Z-bursts taking place within 
the GZK zone of « 50 Mpc. 

This hypothesis was discussed in several pa- 
pers T^-|TJ]. In Ref. pi] , particle spectra were deter- 
mined numerically for case studies which supported the 
Z-burst scenario. The required UHE^ fluxes for differ- 
ent spectral indices were calculated in Ref. |L2|], too. 
The effect of possible lepton asymmetries was studied in 
Ref. |l3| . In Ref. jl5| , the analysis of the Z-burst mech- 
anism was advocated as one of the few possibilities for 
an absolute v mass determination and its potential com- 
pared to others like e. g. the [3 decay endpoint spectrum 
and the i^-less (3f3 decay. 

There is now rather convincing evidence that vs have 
nonzero masses (cf. |l7]]). This evidence comes from 
v oscillation measurements with typical mass splittings 
V 5m 2 ~ 10~ 5 4- 0.4 eV. Neutrinos in this mass range 
are important cosmologically since they represent a non- 
negligible contribution to dark matter (DM) which im- 
poses upper limits on v masses UM- Hydrodynamic 



simulations with massive i^s and including recent ob- 
servational measurements and cosmological constraints 
give (19) J2i m ^ ^ 2 - 4 ' (fiAf/0.17 - 1) eV, if the matter 
content of the universe J7m is assumed to be between 0.2 
and 0.5, as favoured by recent measurements (cf. ]l7||). 

77. Z-burst spectrum and UHECR data. — Our com- 
parison of the Z-burst scenario with the observed UHE- 
CR spectrum is done in four steps. First we determine 
the probability of Z production as a function of the dis- 
tance from earth. In the second step we exploit collider 
experiments to derive the energy distribution of the pro- 
duced protons in the lab system. The third ingredient is 
the propagation of the protons, i. e. the determination 
of their energy loss due to pion and e + e~ production 
through scattering on the CMBR and due to their red- 
shift. The last step is the comparison of the predicted 
and observed spectrum and the extraction of the mass of 
the Ris and the necessary UHE^ flux. 

For a given neutrino type i the probability of Z-bursts 
at some distance r is proportional to the number den- 
sity n Vi (r) of the Ri^s and to the flux F Vi (E Vi , r) of the 
UHE^s at energy E Vi « E T J: S . The density distribution 
of R^s as hot DM follows the total mass distribution; 
however, it is expected to be less clustered. This is the 
reason why we, similarly to Ref. |l4| but in distinction 
to practically all previous authors [10|-^2|, do not fol- 
low the assumption of having a relative overdensity of 
/„ = 10 2 -i- 10 4 in our neighbourhood. For distances be- 
low 100 Mpc we varied the shape of the n Vi (r) distribu- 
tion between the homogeneous case and that of m to t ( r ) , 
the total mass distribution obtained from peculiar veloc- 
ity measurements po|| . Our results are rather insensitive 
to these variations. Their effect is included in our error 
bars. For scales larger than 100 Mpc the R^ density is 
given by the Big Bang cosmology, n Vi = 56- (1+z) 3 cm" 3 . 
In our analysis we go up to distances of redshift z = 2 
(cf. pH), We include uncertainties of the expansion rate 
(see e. g. Sect. 2 of 0). The UHEi/ flux is assumed to 
have the form F Vi {E Vi ,r) = F Vi (E Vi , 0) (1 + z) a , where 
a characterizes the source evolution (see also P,|ll|). In- 
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dependently of the production mechanism, v oscillations 
result in a uniform F Ui mixture for the different types i. 

The Z-burst scenario is based on Z decays. At LEP 
and SLC millions of Z bosons were produced and their 
decay analyzed with extreme high accuracy. 69.89% of 
the Z decays are hadronic and the p + p multiplicity is 
(N p ) = 1.04 ± 0.04 in the hadronic channel 0. The 
neutron multiplicity, which we included in our analysis, 
is 1=3 4% smaller than the proton's |22]. We combined 
existing published and some improved unpublished data 
on the momentum distribution V (x — p pr oton/Pboam) of 
protons in Z decays |p3|j . Due to the large statistics, the 
uncertainties related to Z decay are negligible. 

In the CM system of the Z production the angular 
distribution of the hadrons is determined by the spin 
1/2 of the primary quarks and thus proportional to 
1 + w 2 = 1 + cos 2 9 (here 9 is the angle between the 
incoming neutrinos and the outgoing hadrons (cf. pi])). 
The energy distribution Q(E p ) of the produced protons 
with energy E p is finally obtained after a Lorentz trans- 
formation from the CM system to the lab system, 
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±y - w^/y 2 - (1 - w 2 ){2m p /M z ) 2 



y/y 2 -(l-w 2 )(2m p /M z ) 2 
([-toy ± \/y 2 -(l-w 2 )(2m p /M z ) 2 }/{l - w 2 )^j 



where m p is the p mass and y = 2E p / E v . 

Particles of EG origin and energies above « 4 • 10 19 eV 
lose a large fraction of their energies |lj . This can be de- 
scribed by the function P(r, E p , E), the probability that 
a proton created at a distance r with energy E p arrives 
at Earth above the threshold energy E |2j| . It has been 
calculated for a wide range of parameters in Ref . |^6| , and 
the respective data are available at www.desy.de/~uhecr. 
Note, that the energy attenuation length of 7s are longer 
than that of protons roughly by a factor of 10 at su- 
perhigh energies such as 10 21 eV. The detailed study of 
the boosted Z-decay (data from Ref. Q) results in 7s 
of energy below 10 19 (1 eV/m^) eV, where their atten- 
uation length is much smaller, for strong enough radio 
background. Thus, their contribution to the UHECR 
spectrum is far less relevant than that of the protons. 

The Z-burst contribution to the UHECR spectrum, for 
degenerate v masses (m„ ss m Vi ) , is given by 

/>oc pRo z" 00 

j{E,m„)=I-F z 1 ■ dE p dr de (2) 
Jo Jo Jo 

F Vi {E Vi , r)a(e)n Vz (r) Q(E P ) (-dP(r, E p , E)/dE) , 

i 

where / w 8 • 10 16 m 2 -s-sr is the total exposure (es- 
timated from the highest energy events and the cor- 
responding fluxes), Rq is the distance at z = 2, and 



a(e) is the Z production cross section at CM energy 
e = (2m„ E Ui ) x / 2 , The normalization factor Fz is pro- 
portional to the sum of the v fluxes at CM energy M z . 

We compare the spectrum (|^) with the observed one 
and give the value of m v based on a maximum likelihood 
analysis. In the Z-burst scenario a small Ris mass needs 
large El es in order to produce a Z. Large E™ s results in 
a large Lorentz boost, thus large E p . In this way the 
detected E determines the mass of the Rza 

Our analysis includes the published and the unpub- 
lished (from the www pages of the experiments on 
17/03/01) UHECR data of §-|,|. Due to normaliza- 
tion difficulties we did not use the Yakutsk Q results. 

Since the Z-burst scenario results in a quite small flux 
for lower energies, the "ankle" is used as a lower end for 
the UHECR spectrum: log(-E min /eV) = 18.5. Our re- 
sults are insensitive to the definition of the upper end 
(the flux is extremely small there) for which we choose 
log(£' max /eV) = 26. As usual, we divided each logarith- 
mic unit into ten bins. The integrated flux gives the total 
number of events in a bin. The uncertainties of the mea- 
sured energies are about 30% which is one bin. Using 
a Monte-Carlo method we included this uncertainty in 
the final error estimates. For the degenerate case, the 
predicted number of events in a bin is given by 



N(i) 



dE [A-E-P + Fz-j{E,m u )] 



(3) 



where Ei is the lower bound of the i th energy bin. The 
first term is the usual power law, which describes the 
data well for smaller energies (|]. For this term we will 
study two possibilities. In the first case we assume that 
the power part is produced in our galaxy. Thus no GZK 
effect should be included for it ("halo"). In the second 
- in some sense more realistic - case we assume that 
the protons come from uniformly distributed, EG sources 
and suffer from the GZK cutoff ( "EG" ) . In this case the 
simple power-law-like term will be modified and falls off 
around 4 • 10 19 eV (see later Fig. |l|). The second term 
of the flux in Eq. (||) corresponds to the spectrum of the 
Z-bursts, Eq. (|^). A and Fz are normalization factors. 

The expectation value for the number of events in a 
bin is given by Eq. (||) and it is Poisson distributed. 
To determine the most probable value for m„ we used 
the maximum likelihood method and minimized |27|] the 
X 2 (/3, A, Fz, vn v ) for Poisson distributed data (IT; 

26.0 

X 2 = 2[JV(t)-JV (i)+JV (t)ln(JV (*)/iV(*))], (4) 

i=18.5 

where N a (i) is the total number of observed events in the 
z th bin. In our fitting procedure we have four parameters: 
f3,A, Fz and m„. The minimum of the ^{(3, A, Fz, m v ) 
function is Xmin a ^ ~*Vmin> the most probable value for 
the mass. The 1 <r confidence interval for m„ is given 
by X 2 (P[, A', F' z ,m v ) = Xl 2 nin +1. P' ,A> ,F' Z are defined by 
minimizing x 2 (Pj A, Fz, vtiy) in (3, A and Fz at fixed m„. 
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FIG. 1. The available UHECR data with their error bars 
and the best fits from Z-bursts. Note that there are no events 
above 3 x 10 20 eV (shown by an arrow). The dotted line shows 
the best fit for the "halo" -case. The bump around 4 ■ I0 19 eV 
is due to the Z-burst protons, whereas the almost horizon- 
tal contribution is the first, power- law-like term of Eq. (|^). 
The solid line shows the "extragalactic"-case. The first bump 
at 4 ■ 10 19 eV represents protons produced at high energies 
and accumulated just above the GZK cutoff due to their en- 
ergy losses. The bump at 3 • I0 21 eV is a remnant of the 
Z-burst energy. The dashed line shows the contribution of 
the power-law-like spectrum with the GZK effect included. 
The predicted fall-off for this term around 4 • f 19 eV can be 
observed. The attenuation of the Z-burst component appears 
to be weaker on account of the narrowness of the injected 
proton spectrum and the fact that the observed post-GZK 
protons are produced within the GZK zone. 



Our best fits to the observed data can be seen in Fig. [j], 
for evolution parameter a = 3. The neutrino mass is 

2J5 -a97(i'.89) eV for the "halo"- and 0.26^;u(a22) eV 
for the "EG" -case, respectively. The first numbers are 
the 1 a, the numbers in the brackets are the 2 a errors. 
This gives an absolute lower bound on the mass of the 
heaviest v of 0.06 eV at the 95% CL. Note, that the 
surprisingly small uncertainties are based on the above 
X 2 analysis and dominantly statistical ones. The fits 
are rather good; for 21 non-vanishing bins and 4 fit- 
ted parameters they can be as low as x 2 = 18.6. We 
determined m v for a wide range of cosmological source 
evolution (a — 0^3) and Hubble parameter (Hq = 
(71 ± 7) xj^gg km/sec/Mpc) and observed only a moder- 
ate dependence on them. The results remain within the 
above error bars. For these mass scales the atmospheric 
or solar v experiments suggest practically degenerate v 
masses. This has no influence on our v mass determina- 
tion, but is taken into account in our flux determination. 

We performed a Monte-Carlo analysis studying higher 
statistics. In the near future, Auger |28| will provide a 
ten times higher statistics, which reduces the error bars in 
the neutrino mass to « one third of their present values. 
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FIG. 2. The best fit and the 1 a (68% confidence level) 
region in a scenario with two non-degenerate v masses. 



One of the most attractive patterns for v masses is 
similar to the one of the charged leptons or quarks: the 
masses are hierarchical, thus the mass difference between 
the families is approximately the mass of the heavier par- 
ticle. Using the mass difference of the atmospheric v os- 
cillation for the heaviest mass (l7), one obtains values 
between 0.03 and 0.09 eV. It is an intriguing feature of 
our result that the smaller one of the predicted masses is 
compatible on the ss 1.3 ct level with this scenario. 

Another popular possibility is to have 4 neutrino types. 
Two of them - electron and sterile neutrinos - are sepa- 
rated by the solar v oscillation solution, the other two - 
muon and tau - by the atmospheric v oscillation solution, 
whereas the mass difference between the two groups is of 
the order of 1 eV. We studied this possibility, too. On our 
mass scales and resolution the electron and sterile neu- 
trinos are practically degenerate with mass m\ and the 
muon and tau neutrinos are also degenerate with mass 
TO2 . The best fit and the 1 a region in the m\ — mi plane 
is shown in Fig. | for the "EG" -case. Since this two-mass 
scenario has much less constraints the allowed region for 
the masses is larger than in the one-mass scenario. 

III. Necessary UHEv flux. — The necessary UHE^ flux 
at El° s can be obtained via Eqs. (||) and (g) from our 
fits. We have summarized them in Fig. ||, together with 
some existing upper limits and projected sensitivities of 
present, near future and future observational projects. 
The necessary v flux appears to be well below present 
upper limits and is within the expected sensitivity of 
AMANDA, Auger, and OWL. Clearly, our fluxes are 
higher than the ones found in Ref. [jllj based on local 
overdensities f v . However, since we also have a back- 
ground the normalization of the Z-burst component is 
different and correspondingly our fluxes are somewhat 
less than a factor of /„ higher. An important constraint 
for all top-down scenarios M is the EGRET observation 
of a diffuse 7 background |34j . As a cross check, we calcu- 
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[2] M. Takeda e t al, Phys. Rev. Lett. 81, 1163 (1998); astro- 
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FIG. 3. Neutrino fluxes, F = § Y? i= i( F ->i+ F *i)> required 
by the Z-burst hypothesis for the "halo" and the "extragalac- 
tic" case, for evolution parameter a — 3. The horizontal 
errors indicate the 1 a uncertainty of the mass determination 
and the vertical errors include also the uncertainty of the Hub- 
ble expansion rate. The dependence on a is just of the order 
of the thickness of the lines. Also shown are upper limits from 
Fly's Eye ^| and the Goldstone lunar ultrahigh energy neu- 
trino experiment GLUE [|30tl , as well as projected sensitivities 
of AMANDA Hi, Auger |y,|| and OWL 



lated the total energy in 7s from Z-bursts. We assumed 
that all energy ends up between 30 MeV and 100 GeV. 
Our 7 flux is somewhat smaller than that of EGRET. 

IV. Conclusions. — We compared the predicted spec- 
trum of the Z-burst hypothesis with the observed UHE- 
CR spectrum. We should emphasize, that only a real- 
istic overdensity of Ri/s was used. We determined the 
mass of the heaviest R^: m v — 2.75 jlp g® eV for halo and 
0.26^q'^4 eV for EG scenarios. The second mass, with 
a lower bound of 0.06 eV on the 95% CL, is compatible 
with a hierarchical v mass scenario with the largest mass 
suggested by the atmospheric v oscillation. The neces- 
sary UHEz/ flux should be detected in the near future. 
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